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Abstract 

We simulate a network of routers and M network users making 
concurrent point-to-point connections by buying and selling router capac- 
ity from each other. The resources need to be acquired in complete sets, 
but there is only one spot market for each router. In order to describe 
the internal dynamics of the market, we model the observed prices by A'^- 
dimensional Ito-processes. Modeling using stochastic processes is novel in 
this context of describing interactions between end-users in a system with 
shared resources, and allows a standard set of mathematical tools to be 
applied. The derived models can also be used to price contingent claims 
on network capacity and thus to price complex network services such as 
quality of service levels, multicast, etc. 

1 Introduction 

To be able to provide guaranteed quality of service, QoS, in a network, a user 
needs to be able to reserve capacity, or 'bandwidth', in congested routers. The 
reservation scheme should be efficient in the sense that one reservation should 
not unnecessarily block other reservations, and it should not require extensive 
negotiation. One scheme that fits these requirements is to trade router capacity 
in spot markets. The assumption is that someone reserves capacity for a con- 
nection by buying the capacity in the routers along the cheapest path between 
the source and destination node. When the capacity is no longer needed, it is 
sold to someone else. Increased demand increases price, so alternative paths, 
if they exist, may become competitive, and the users will tend to move their 
bandwidth usage away from congested routers. 

We ultimately want to be able to price contingent claims on resources, for- 
malized as options or futures. The working hypothesis is that adding a suitable 
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set of such claims will improve the efficiency of the resource allocation, as prices 
will better reflect all available information. Indeed, the purpose of contingent 
claims on future value, such as options, futures, or more generally, derivative 
securities, is to construct a market for trading the expected future value. For 
the problem studied here, this means trading in expected future demand and 
supply. One way to price derivatives is to use statistical modeling of the price 
dynamics, since under suitable assumptions, which we will cover in a forth- 
coming separate contribution, d erivative prices are func tions of current market 
prices and the statistical model] Biack73| | Avellaneda99 1 . This paper addresses 
the necessary preHminary issue of how to estimate parameters in two stochastic 
models from observed market prices, and also presents measures of the efficiency 
of the market resource allocation scheme. These questions are also of indepen- 
dent interest, as the efficiency is to be used to compare different market system 
with each other, and with other allocation schemes. 

Using artificial markets for resource allocation in distributed sys tems date s 

back to th e mid 80ies , ranging over m arkets for storage cap acity [Kurose8!: |, 

CPU time | Ferguson88 1 | Waldspurger92| , and network capacity | Kurose85 1 |^airamesh95| . 
The emphasis has been on evaluating the efficiency of the resource allocation, 
rather than understanding the resulting price dynamics. More recent work 
stresses the agent aspect, i.e. that the trading parties are locally optimizing 
entities [FaratinOC]. Combinatorial markets, i.e. trading of bundles of distinct 
resources is yet a relatively new area. Somewhat related is the area of combina- 
torial auctions | Rassenti82| |Rothkopf98 1 | ^andholm9S |. The use of derivatives 
for network admission control has been used in | Lazar98| . 

Previous bandwidth market models usually only include a primary market, 
in which end-users can buy and sell capacity only from the router owner. In the 
presented model end-users trade amongst each other, i.e. the capacity is traded 
on secondary markets. 

This paper is organized as follows: in section ^ we present the model, of 
which a central ingredient is a detailed mechanism of price formation, with 
very low overhead in user-user communication. In section H] we model math- 
ematically the resulting price processes. The main tools used are stochastic 
differential equations (Ito-processes), their associated Fokker-Planck equations, 
and the stationary distribution of those. In section ^ we sum up an discuss our 
results. 



2 Method 



2.1 Market model 



We simulate a network consisting of N inter-connected routers and M network 
users concurrently making reservations of router capacity for point-to-point con- 
nections. There is one spot market per router, in which users trade router 
capacity. 



We use Farmer's non-equilibrium market dynamics | FarmerOC | as a prescrip- 
tion of how prices change due to trading (see also sec. Farmer's model is 
based on the assumption of a market maker that guarantees liquidity at all 
times, and that buying causes prices to increase, while seUing causes prices to 
decrease. Farmer assumes that the price dynamics is such that it is impossible 
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Figure 1: The network with 10 routers. 



to move the market by performing a sequence of trades, where the net traded 
volume is zero, and that the relative price changes are independent of the cur- 
rent value of the price. From these assumptions. Farmer derives a formula that 
the price per unit in a transaction of w units is S{S,uj) = Se'^/^, where S is the 
unit price in the previous trade, and A the market depth or Hquidity, i.e. the 
rate at which the price is changed by trading. For a derivation, see Appendix A. 
This model is useful since we do not have to simulate details of the order-book 
in each market. Instead, we can calculate the price change caused by trading 
directly from the last price and the size of the trade. 

2.2 Simulation setup 

At each time interval, m new demands are generated. A demand is a 7-tuple, 

d=< id, uid, src, dst, cap, dur, max > 

The user identities are chosen independently, so one user may receive zero or 
one or several new demands. If the current time is t, demand id specifies that 
the user uid demands cap units in each node on a path from src to dst, starting 
at time t and ending at time t + dur if it costs less than max to obtain the 
resources. If there are several paths, a choice will be made, see below. During 
the simulation a user reserves capacity in a router by buying that capacity, and 
sells excess capacity that is no longer in need. 

User i owns r,j units of j. Initially, none of the simulated users are assigned 
any resources or money, i.e. rij = 0, cashi = for all i and j, nor do they 
have any resource demand, ivij = 0. When a user manages to satisfy a demand, 
its capital is increased by the amount max, and when resources are bought 
and sold, the capital is decreased and increased, respectively, by the cost of the 
resource. The simulation is run in L time steps from time to T with time 
increments At. At each time step where the current time is t: 

• Generate m new demands. A demand is specified by a unique demand 
number id, a user uid randomly drawn from the set of M users, a source 
node src and a destination node dst, both randomly drawn from the set 
of N nodes, and the required capacity per node cap, which is randomly 
drawn as the ceiling of e''^^ where ^ is a uniform random variable between 



3 



and 1. If some demand is more important than others, a user is willing 
to pay more for that resource. In this simulation, max is a linear function 
of the required capacity, i.e. the maximum total cost max = Cunit cap, 
where C'unit is a simulation parameter. The duration (number of time 
steps) dur is randomly drawn integer between 1 and D. 

• Calculate ujij , the net change of resource j of user i — uid in the following 
way: 

- For each new demand d, user i looks at the last known transaction 
prices Sj{t) and decides to buy cap resources along the least cost 
path. Since the user will not know the actual cost of a buying and 
later selling resources along a path, its decision on whether to buy 
resources or not, is determined by a parameter Cmax- The user de- 
cides to buy the resources if and only if the estimated total cost to 
buy the resources is less than Cmaxmax. The resulting demand cap 
is then added to the LOij for all the resources j on the least cost path, 
and the amount max is added to cashi. 

- For each satisfied demand d terminating at t, decide to sell the re- 
sources that were allocated (if any) to the demand, i.e. subtract the 
resulting supply of liberated capacity on the least cost path if the 
demand was satisfied and router capacity was bought. 

• All the demands on capacity in single routers (i.e. ujij) are effectuated 
one by one in random order, and for each trade prices are updated to 
Sjc'^'-'^/'^' . The cashi is decreased by (^ijSj and the owned amount of 
resource is increased by uJij. According to the price formation formula, 
the trades are made at prices which depend on the actual order. The prices 
payed by the users are thus not the same as those used for determining 
the least cost path. The next known price is the one that pertains after 
all the trades, and is independent of the order. 

• When all trading is done for this time-step, log the last transaction prices, 
S{t), the number of satisfied demands, and repeat. 

2.3 Simulation Parameters 

The simulation was run for L = 1000 time-steps {At — 0.01) using the network 
in (fig. 1^). There were 10 routers and 10 users, so iV = 10 and M — 10. The 
liquidity in all Farmer markets was chosen to be = 10. The maximal cost per 
route users were willing to accept was determined by Cunit = 100 and Cmax — 1 
for all users. Initial prices were set to Si{0) = 10 in all markets. Every time 
step, m = 10 new demands were generated. The duration dur was uniformly 
distributed between 1 and 10, i.e. D — 10. The call duration sets the time scale 
in the simulation, as we will see below. The required amount of capacity was 
determined hy K = 2 (see above). 

3 Results 
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Figure 2: Plot of the unit price of capacity in router 1 from time 0.01 to 10. 
At=OM. 

3.1 Statistical modeling of the price process 

Inhomogeneous drift, additive noise, constant coefficients The price 
as shown in (fig. ||) does not appear to be drifting freely. Instead it appears 
to return towards the same area. Assuming the price process is an Ornstein- 
Uhlenbeck process, the dynamics for the price of router i would be 

dS.it) = a,{^i, - S,it))dt + a,dW,{t) (1) 

where Wi{t) is a Wiener process and the correlation between two processes i and 
j is Corr[dWj{t), dWj{t)] — pij. Recall that a Wiener process has independent 
normal distributed increments with mean and variance t — s, i.e. Wi{t) — 
W{si) ~ A^[0, y/t — s] where t > s. For improved readability, we omit the 
indices in Si{t), etc., when they are irrelevant for the understanding. 

In (eq. |l|) the drift term (the dt term) detracts when S is bigger than /z, 
and adds when S is less than /i. The amount of the increase is determined by 
a. The diffusion (the dW term) is independent of S. 

When the simulation has run for sufficiently long time, the price Si{t) be- 
comes independent of the starting state of the system, and reaches a stationary 
distribution P{t) of S{t). Using the Fokker-Planck equation (see Appendix A) 
gives 

P{s) = Cae ^^"/v^' (2) 

where Co = {'na'^ /a)^^ , which is the density function for a normal distribution 
N\u,,-?=\. We note first that the normal distribution is non-zero in all of 

v2a 

(— oo, oo), meaning that S{t) can take on negative values, something that is not 
possible in a market with Farmer's dynamics. Second, the normal distribution 
has many good properties, such as that a weighted sum of normal distributed 
variables is also normal distributed. If prices are far from zero, this dynamics 
may therefore be a convenient approximation to the true distribution. Third, 
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the stationary distribution only depends on the ratio ^ and can therefore not 
distinguish between separate variations in and a, which has to be done by 
other means. 

The observations S(i), i e [1, ■■■,L] of the process S in its stationary state 
are regularly spaced with distance At. Since E[S] — ji, we estimate ^ with 



4 = 1 

We estimate a after noting that E[{dSY] = a^dt + 0{dtdW). Therefore, for 
small At, 
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Since yar[S'] = the unbiased estimate is 

L 



2d L-1 

so the estimate of a is 



E[S'] -E[S]' 



2At 



{J:tlS{^r-{EllSi^))'■ 



To verify that the observed estimate is indeed correct, we can plot the his- 
togram of the observed data together with the estimated density function P{s), 
see (fig ^ . This estimation coincides with the maximum likelihood estimate of 
the variables. A least square fit to the values in the histogram gives a better- 

2 

looking curve. However, this fit only gives ^ and the ratio as noted above. 

Having estimates of a, /i and a we are able to estimate the correlation pij 
between the random sources for two price processes Si and Sj by solving dWi 
and dWj in (eq. |]) and using the estimates 

n ^ ^*(^ + 1) - ^«(^) - - ^4(fc))At 
dW^{k) = i 



so that the estimate becomes 



Gov[dW^,dWJ] 
YaT[dw,]'VaT[dWj] 
T^iJ:tldm{k)dW,{k) 



(At)2 

where we have used = and yar[o?VF] = dt. 
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Figure 3: Histogram over the observed prices, the estimated normal distributed 
density function N[^,(j/y/2a] (the red dashed Une), and a least square error fit 
to the 15 columns in the histogram (the green dotted line). 



Inhomogeneous drift, multiplicative noise, constant coefficients The 



dynamics in (sec. |3.l[) can generate negative prices, something which is should 
not be possible in a well-functioning market. By asserting a mean-reverting 
dynamics with multiplicative noise we get a process that is strictly positive. 
Assume the dynamics for the price of router i to be 

dS,{t) = a,{^l, ~ S,{t))dt + cjA{t)dW,{t) (3) 

where Wi{t) is a Wiener process and the correlation between two processes i 
and j is Corr[dWj{t), dWj{t)\ = pij. As before, we we omit the indices in Si{t), 
etc., for readability. 

The stationary distribution P{s) of S is 

= (4) 



r(7) 

where 7 = and T{z) is the gamma function. In (eq. H) s takes only positive 
values, which is consistent, since the dynamics in (eq. does not move an s 
from positive to negative. 

The first moment of this distribution, E[S] is /i (see Appendix A). As for 
the Ornstein-Uhlenbeck process above, we estimate with 

A = ^[5] = -^^(z) 

We cannot easily estimate a from E[{dS)'^] since S is in the drift term, but note 
that the process X{t) — log{S{t)) in (eq. @) has additive noise. Therefore, we 



7 



get for small At, 



1 

At 



1 1 ■^^"^ 



S{i + 1) 



To estimate a, we try two approaches. First we use the conditional expectation 
E[S{t + T)\S{t)] for S in the stationary state of the system. Taking the partial 
derivative w.r.t r gives a first order ODE with the solution 



E[S{t + T)\Sit)] 



showing that the expected value approaches exponentially with r and a de- 
termines the speed of the return. Rearrange to keep e~"'^ (which is independent 
of S{t)) on one side, take the logarithm and take the expected value of both 
sides. Let r = fcAt. We now have an estimate for a. 



exp{^a kAt) 



1 



L-k 



L-k 

E 



S(i + k) — pL 



S{i)-[i 

As we can see from the left hand side of the equation above, plotting the right 
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Figure 4: Left fig.: the estimate y = e~°-^^^ derived from the conditional expec- 
tation of S'(t-l-r) for the first 20 time-steps in blue, and the plot of y = e""'^^* in 
green. The estimate may be less than zero due to the noise from the simulation. 
Right fig.: the covariance between S(t -I- fcAt) and S(t\ The auto-correlation 
deviates some from an exponential, possibly because the price dynamics has 
higher order than assumed. It would explain why the estimate in the left plot 
deviates from an exponential. Simulation time, 100 000 time steps, — 10, and 
D = 100. 



hand side as a function of fc should result in a straight line if a is constant. 
However, as we can see in (fig. ^, the fine is straight up to some fc that depends 
on the simulation parameters and error. The estimation is very sensitive to 
errors in /i, especially in the denominator if S is near /i, which results in a 
fiattened out jagged line, since e~°''"^*' is less than the simulation error for 
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Figure 5: Histogram over the observed prices and the estimated gamma distri- 
bution. ^ and a are estimated using the moment method, and after that, a is 
estimated with a least square error fit. 



larger k. It eventually flattens out which means that a appears to decrease for 
larger t. 

Another way to estimate a is to estimate fi and a as above, and then fit the 
observed distribution to the distribution (eq. ^ using the least square method. 
See (fig. H) for a plot of the model fit. 

Having estimates of a, ^ and a we are able to estimate the correlation 
between the random sources for two price processes in the same way as above, 
by solving dWi and dWj in (eq. ^) and using the estimates 



AM^,(fc) = 



S,{k + 1) - S,{k) - - S,{k))M 



atSi{k) 



so that the estimate becomes 



zhEti dW^{k)dW,{k) 

(At)2 



where again we have used i?[c?VF] — and yar[(iW^] — dt. 

Plotting the histograms of ISW shows us if the model is good. If that is the 
case, the AM^ should behave like samples of a Wiener process, i.e. be normal 
distributed. As can be seen in (fig. ||), the model (eq. ^) fits well if the market 
liquidity is high. However, if market liquidity is low then the assumed price 
model does not provide a good fit. 



Price Correlations The resource prices in a network of resources depend 
on the prices of other resources, since they are traded in groups. Using the 
parameter estimates of the price dynamics above, we get the correlation matrix 
in (tab. |l|) for a simulation with high liquidity markets (Ai = 100). 

We find that prices generally are positively correlated. Comparing with the 
network graph (fig. |l|), one can see that prices of neighboring nodes often are 
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Figure 6: The left figure shows the histogram (in red) over ISW in the simu- 
lation described in (sec j2.3| ) with \i = 10, and a fitted density function for a 
normal distribution (in blue) . The right figure shows the histogram of another 
simulation with the same parameters except for the market Hquidity which is 
A, = 100. 

strongly positively correlated with an average correlation coefficient of around 
0.4, compared to 0.03 for nodes that are not connected. The most significant 
exception is nodes 3 and 7. Looking at the network graph it is clear that no 
least cost path can contain both 3 and 7, except for the path from 3 to 7. 

3.2 Efficiency of the market based resource allocation 

The efficiency of a market based resource allocation scheme depends on how well 
prices reflect the available information about resource demands. Two sources 
for bad performance is that price quotes are outdated, or that they do not 
reflect knowledge about the price behavior, e.g. periodic price fluctuations, or 
future prices. To measure efficiency, we can measure the resource utiHzation to 
compare which markets are able to capture most information. 

Communication costs The number of messages sent in the negotiation phase 
of this kind of system is negligible, since all users operate on old price quotes 
and place bids at market, i.e. accept the price whatever it may be. They do not 
update their price quotes for every bid. Therefore, less than M messages per 
trade come in to a user (the potential quote update from each market). One 
message per trade (the bid) go out from each user to the markets that contain 
resources that the user have chosen. No messages need to be communicated 
between the end-users. 

Successful connection ratio In the Farmer market model, the less liquid the 
markets are, the less valid are the price quotes. To the left in (flg. 0) the ratio 
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Table 1: Correlation coefficients for AW for a simulation with Xi — 100, after 
5000 time steps. Connected routers are in bold face. 



of successful connections is plotted as a function of t he market liquidity for 
simulations with the parameters described in (sec. O). The graphs corresponds 
to different values of the decision parameter C'max ■ Low liquidity causes large 
price fluctuations, making the prices higher than the Hmit CmaxTnax (see sec. 
2.2) which inhibits many connections. Note that a connection is considered 
successful even if the net cost (after releasing the resources) is higher than max. 



Net profit To the right in (fig. |^ we plot the average profit as a function 
of the liquidity for a number of values of Cmax- Large values of Cmax causes 
the users to buy resources when they are expensive. If the liquidity is high, 
the users sell resources at approximately the same price, but with low liquidity, 
prices will move significantly downwards when the resources are sold, causing a 
net loss to the user. 

The particular way the max cost for a connection is determined of course 
very much determines which of the different kinds of traffic that is promoted in 
the network. Different schemes could for instance promote short or long paths, 
high or low capacity connections, etc. 



Average load The average load, or reserved capacity in a router is with the 
Farmer dynamics a direct function of the market price, Ljj{t) — Alog^^. If two 
simulations that allocate the same connections differ in load, the higher load 
depends on inefficient routing that does not choose efficient path. If the simu- 
lations differ in allocated connections, the comparison is less straight forward, 
and depends on the choice of metric above. We intend to return to a longer 
discussion on efficiency, in particular for load balancing, in a future contribution. 



4 Discussion 

In our simulation. Farmer's market dynamics can be interpreted either as pre- 
scriptions of how a rational market maker should modify the market prices, or as 
a model of the aggregate behavior of market price changes during some period 
of time. 



11 




100 



Figure 7: To the left, average ratio of successful connections as a function of the 
market liquidity. The graphs from bottom to top have Cmax equal to 1, 4, 16 
and 64. To the right, average net profit as a function of market liquidity. The 
graphs from bottom to top have Cmax equal to 64, 16, and 4, showing that a 
high Cmax causes losses on average in a low liquid market. 

With the former interpretation, the Farmer dynamics can be used to im- 
plement a market maker program that brokers trading between end-users. The 
dynamics is derived from the assumption that it is impossible to change prices 
by 'trade in circles'. This assumption is a necessary condition for any mar- 
ket maker strategy, since otherwise anyone can exploit the market and gain an 
unlimited amount of money from the market maker. 

The latter interpretation can be used when we want to simulate a part of 
a market with many concurrent trades. Since trading causes prices to change 
means that it is impossible to have updated price information at the time they 
place their bid. It is only possible to have completely updated price info if the 
trading is synchronous, something that severely reduces the number of bids a 
market can handle per time unit if communication delays are taken into account. 

The market Hquidity (or market depth) parameter Ai determines the speed 
at which prices change. When the dynamics interpreted as prescriptions for a 
market maker, it is up to the market maker to adjust in order to reduce the 
risk of running out of resources. If we model an existing bandwidth market, Ai 
is an observable parameter which must be determined from the distribution of 
price jumps in relation to traded volume. The more resources in relation to the 
order size, the larger A^. A bandwidth market that is trading the capacity of 
a high capacity router can, everything else being equal, be expected to have a 
higher Ai than a similar market for a congested or inefficient router. 

In the mean reverting processes investigated here, a determines the speed 
with which the process returns towards its statistical average /i. The speed of 
return depends on the characteristic time length of the system, determined by 
the simulation parameter Z?, the call length. With a large D, the prices will 
be correlated with previous prices over a longer time period, which results in a 
small a. With a small D, the price effect of previous calls will soon be forgotten. 
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resulting in a high a. 

The traditional statistical models used for telephony have been found to in- 
adequately model data network traffic. Data communication has been found 
to show a very bursty or fractal behavior as one communication event often 
generates a burst of more communication to other parts of the network. Data 
communication is generally short lived and often with strong latency bounds, 
due to the increasing use of computer networks for interactive communication. 
For short lived communication with low transfered amount, traditional switched 
best effort networks will probably continue to provide a very efficient solution 
for a long time. However, with increasing demand of streaming real-time data 
such as high quality video-telephony, video-on-demand, etc., it is necessary to 
be able to reserve capacity. The alternatives are large buffers, which has bad 
latency performance, or migrating data (intelligent replication), which can re- 
duce load for one-to-many communication. Neither is suited for point-to-point 
communication service guarantees. 

In the 'bandwidth market' presented in this paper, end users buy the required 
resources themselves. The demand for updated price quotes results generates 
additional network traffic. In an extended model, we could allow risk neutral 
middle men to sell options on the resources. Since updated quotes will then only 
be required by the (few) middle men doing the actual trading, the additional 
traffic would be greatly reduced. We intend to return to these topics in future 
work. 
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Appendix A 
Stationary Distribution 

S{t) is an Ito process with the dynamics dS{t) = a{t, S{t))dt + h{t, S{t))dW{t), 
and P{s, T; sq, t) is the contingent probability distribution of S at time T given 
S{t) = So- P obeys the Fokker-Planck equation a.k.a the forward Kolmogorov 
equation, 

- OtP - ds (a{T, s) P) + a,, (^^^P) = (5) 
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An important class of stochastic processes are those that are stationary. For 
those, if T — t is sufficiently large, the contingent distribution P no longer 
depends on sq, and can be substituted with the stationary distribution P{s), 
which does not depend on so, t or T. Two ways in which a stochastic process 
can fail to be stationary is if the coefficients, e.g. a and b, are explicitly time- 
dependent, or if there is no distribution P{s). The latter happens for instance 
in ordinary random diffusion, in which the probability gradually spreads out 
without reaching a limit. 

Assuming that a stationary distribution has been reached, the time deriva- 
tive drops out of (eq. ||), and the equation can be integrated once to be 



2a{s) - 2b{s)dsb{s) 



which is an ODE. One further integration of (eq. ||) gives 

PGs) = Ceio W^"^ (7) 

where C is a normalization constant determined by / P{s)ds — 1. 

The Additive Noise Process 

Stationary density function, mean and variance Assume the price dy- 
namics 

dS{t) = a{p - Sit))dt + adW{t) 
Then the conditional stationary probability distribution of is, using (eq. 0), 

P{s) = CqB Jo ^'^ 



= Coe ^^-/-y^' 

where Co — {na'^ /a)~i . P{s) can be identified as the density function for a 
normal distribution with mean /i and variance a"^ /2a. 

The Multiplicative Noise Process 
Stationary density function Assume the price dynamics 
dS{t) = a{fi - S{t))dt + aSdW{t) 

Using (eq. 0) 

= log S -log SQ 

P(s) = Cie -'-o " ( — ) 

Vers/ 



C2 \S, 

An alternative derivation is to let X(t) = log5'(i). Ito's lemma gives that 

2 

dX = (^afxe^^ - (a -I- y -I- adW (8) 
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Then the stationary probabiHty distribution Q{x) of X is , using (eq. 



) du 



1 



= Cie"'''^"^ "'+(7+1)2: 

where Ci is a normalization constant, and 7=2^ for readabiUty. The stationary 
density function P{s) for S is found by recaUing X — log(S') and 



P{s)ds = g(log(s))d(log(s)) 



The constant C2 is determined by normaHzation and identifying the integral 
as a gamma function, T{z + 1) = t^e~'^dt = zT{z). After the substitution 
t = we get 



/ 1 S7+2 1 



e 



/ 1 \7+l 



00 



^7^^^ Jo 

(— -r(7) 



stationary mean Use the same substitution, t — to obtain 



poo 

E[S] = / sP{s)ds 
Jq 



7Ai 62 V7U / 



00 7Ai (-2 ^7M 
C2^7Ai^ Jo 

=r(7) 

= M 

Conditional expected value Let f{T\S{t)) be the conditional expectation 
of S'(t + T) given S{t). 



fir\S{t)) EE E[S{t + r)\S{t)] 

rt+T 

= Sit)+E[ dSit')\Sit)] 

ft + T 

= S{t)+afiT-aj E[S{t')\S{t)]dt' 

Take the partial derivative with respect to r 

drf{T\S{t)) = afi-aE[S(T)\S(t)] 
= afi-af{T\{S{t)) 
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Multiplying with e""^ and collecting the terms gives 

a,(e"V(T|(S(t))) =e"-aM 

Integrate both sides over t from to t and multiplication with e^""^ gives 

/(r|5(t)) = e-"^(5W-/i)+/x 

where we have used that f{0\S{t)) = S{t). Therefore 

^ f{T\S{t))-n 
Sit) - A. 

The left side does not depend on S(t). Taking the expected value of both sides 
= E 



E[S{t + T)\S{t)]-^i 



f 



S{t) - M 
E[S{t + T)\S{t) =x]- II 



X — jl 



P{x)dx 



where P is the stationary distribution of S, and we have substituted a time 
average with an ensemble average, which is true under suitable assumptions. 

Estimation of parameter a We can estimate the integral from the observed 
data using the estimates 



E[S{t + T)\S{t)=x\ = 



t 



\{x\S{t) = x}\ 



where 5{x) is the Dirac delta function and \A\ is the cardinality of the set A. 
The estimation is 



Jo 



X — n 



= f{j^T.S{x-S{t)))E[^^^^^\S{t)=x]dx 

= -E 



S(t + T)-il 

S{t)-fi 
so we can estimate a with 



which should be a constant function. 
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Farmer's market dynamics 

Let S be the current price, to the net demand and S the price at which the 
demand can be met. We seek a functional relationship of the type S = S{S,uj), 
where S depends continuously on S and oj. Assume the price to be positive 
and bounded, S to be an increasing function of ui, and that prices are only 
changed through trading, S{S, 0) = S. Assume furthermore that one cannot 
make money by trading in circles 

S{S{S{S, uji),LU2),-{i^i + '^2)) - S (9) 

and the relative price change is independent of the absolute price, 

S{S, uj) 



S 



H^) (10) 



Because of ^ we see that S{S{S,uj),—lj) = S, so the inverse is S~^{S,uj) = 
S{S, — w), where is the inverse function of S. Applying 5'~^(-, —{uji + 0J2)) 
on (H) to see that 

S{S{S,u;^),u;2)^S{S,u;,+LU2) (11) 



L0[) together with ( p. l|) implies that </)(wi)(/)(a;2) = (f'i^i + ^2) which implies 
) = e^/^ for some constant A. Therefore 

S{S,uj) ^ Se'^/^ (12) 
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